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More than twice as much carbón is held in soils as in vegetation or 
the atmosphere 1 , and changes in soil carbón content can have a 
large effect on the global carbón budget. The possibility that 
climate change is being reinforced by increased carbón dioxide 
emissions from soils owing to rising temperature is the subject of a 
continuing debate 2-9 . But evidence for the suggested feedback 
mechanism has to date come solely from small-scale laboratory 
and field experiments and modelling studies 2-9 . Here we use data 
from the National Soil Inventory of England and Wales obtained 
between 1978 and 2003 to show that carbón was lost from soils 
across England and Wales over the survey period at a mean rate of 
0.6% yr -1 (relative to the existing soil carbón content). We find 
that the relative rate of carbón loss increased with soil carbón 
content and was more than 2% yr 1 in soils with carbón contents 
greater than 100 gkg -1 . The relationship between rate of carbón 
loss and carbón content is irrespective of land use, suggesting a 
link to climate change. Our findings indícate that losses of soil 
carbón in England and Wales—and by inference in other tempér¬ 
ate regions—are likely to have been offsetting absorption of 
carbón by terrestrial sinks. 

The National Soil Inventory was made to obtain an unbiased 
estímate of the distribution of the soils of England and Wales and of 
the chemistry of the topsoil (0-15 cm depth) 10 . Samples were 
collected and soil profiles described at the intersections of an 
orthogonal 5-km grid over the whole area (Methods). This yielded 
about 6,000 sites, of which 5,662 could be sampled for soil. Figure la 


shows the distribution of soil organic carbón contents across England 
and Wales measured in the original sampling (1978-83). Sufficient 
subsets of the sites were resampled at intervals from 12 to 25 yr after 
the original sampling to be able to detect changes in carbón content 
with 95% confidence (Methods). This was done in three phases: in 
1994-95 for arable and rotational grassland sites (853 of the original 
2,578 sites), in 1995-96 for managed permanent grassland sites (771 
of the original 1,579), and in 2003 for non-agricultural sites (bogs, 
scrub, rough grazing, woodland, and so on; 555 of the original 1,505). 
Roughly 40% of the original sites were resampled. This is the only soil 
inventory on such a scale anywhere in the world to have been 
resampled. To allow for the varying time interval between samplings, 
annual rates of change in carbón were calculated for each site by 
assuming that the process of change was linear over the sampling 
interval. An analysis of known rates of change in soil carbón under 
different conditions showed this to be reasonable. 

Figure 2 summarizes the results grouped by soil type and land use. 
Some differences between soils and land uses are apparent: for 
example, peat soils lost carbón an order of magnitude faster than 
brown soils and man-made soils, and bogs and upland grass lost 
carbón an order of magnitude faster than lowland heath, which 
appears to have gained carbón on average. But we found no 
statistically significant relations between rate of change and land 
use, rainfall class or soil textural class, whether for the data as a whole 
or for outlying data. However, we found a significant negative linear 
correlation between rate of change and original organic carbón 



b Rate of change (g kg -1 yr -1 ) 



Figure 1 Changes in soil 
organic carbón contents 
across England and Wales 
between 1978 and 2003. 

a, Carbón contents in the 
original samplings, and b, rates 
of change calculated from the 
changes over the different 
sampling intervals. Valúes at 
sites that were not resampled 
were calculated from their 
original organic carbón 
contents using equation (1). 
The changes were negative in 
all but 8% of the sites. 
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Figure 2 | Rates of change in soil organic carbón content, grouped by soil 
type and land use. a, Soil type 24 grouping; b, land use grouping. Circles 
indícate mean valúes; error bars indícate 95% confidence intervals. 


content (C org ); that is, the rateofloss increased with C org (Fig. 3). The 
relation applied for the data as a whole and for the three main land 
use groupings separately, though with different slopes and intercepts. 
We took a random subsample of 1,000 observations from the data as 
a whole, and obtained the following relation by linear regression 
specifying an exponential variogram for the residuals (Methods): 

Rate of change in C org = 0.6 — 0.0187 X original C org (1) 

The standard error (s.e.) of the intercept is 0.148 and the s.e. of the 
slope is 0.00081; the rate of change in C org is in units of g kg 1 yr , 
and C org is in g kg 1 . The residuals from this regression applied to the 
whole data set show some regional features, notably a tendency to 
overestimate the rate of loss in uplands in the northwest and South¬ 
west of England. But only 8% of the unexplained variation is spatially 
structured, so a more sophisticated statistical model is unjustified. 
Figure Ib shows the distribution of rates of change across England 
and Wales for all the sites, with valúes for the sites that were not 
resampled obtained using equation (1). 

From the data in Fig. 3, the relative rates of change (rate of change/ 
mean C org over sampling interval, in units of % yr 4 ) were 1.43,0.14, 
— 0.69, —1.84, —2.71, —3.01 and —2.35 for original C org ranges (in 
gkg -1 ) 0-20, 20-30, 30-50, 50-100, 100-200, 200-300 and >300, 
respectively. Henee the relative rates of loss also tended to increase 
with C org , although not above 300 g kg 1 . Soils with C org < 50 g kg 1 
did not lose significant amounts of carbón (that is, not detectable 
over 10 yr), and those with C org < 20gkg 1 appear to have gained 
it. However, for soils with C org > 100 gkg , relative rates of loss 
over the survey period were >2%yr 1 . Given that the bulk of the 


UK’s carbón stocks are in organic soils 11 , this result gives cause for 
concern. 

Table 1 gives estimates of the total changes in carbón in the upper 
15 cm of soil for the whole of England and Wales, with rates for the 
sites that were not resampled predicted using equation (1). The total 
amount of carbón in the upper 15 cm at the time of the original 
sampling is estimated to be 864 Tg, and the total rate of change in this 
depth is — 4.44Tgyr *. From the distribution of soil types in the 
1:250,000 National Soil Map interpolated on a 1-km grid, Bradley 
et al." estímate the total carbón stock of the soils of the UK to be 
2,542Tg over 0-30cm depth (l,015Tg in England, 194Tg in Wales, 
1,161 Tg in Scotland, 172 Tg in Northern Ireland). If the rate of 
change found here also applies to the soils of Scotland (but note that 
most soils in Scotland have large organic carbón contents 11 , and 
therefore the true rate of loss is likely to be larger), and if the rate of 
change over 0-30 cm depth is the same as over 0-15 cm, then pro rata 
the total rate of loss across the UK is 13 Tgyr For comparison, the 
UK’s current industrial C0 2 emission is about 150Tgyr 

We cannot with the existing data say where the missing carbón has 
gone. Some proportion will have been emitted to the atmosphere as 
C0 2 and some will have been leached to deeper soil layers and to 
drainage waters and beyond. The latter would be consistent with the 
observed increases in dissolved organic carbón in surface waters 
across much of England and Wales over the past 40 yr (refs 12, 13). 

Soil carbón contents depend on rates of addition from plant 
growth versus rates of removal in decomposition, leaching and 
other soil processes, and each of these is sensitive to changes in 
land use, climate and other variables 14-17 . Various changes in land use 
will have contributed to carbón losses from soils across England and 
Wales over the survey period, both under agricultural uses (drainage 
schemes, post-war grassland conversión, increased stocking rates) 
and non-agricultural uses (afforestation on wet soils, increased 
erosión, increased burning of upland vegetation). However, we do 
not have sufficient data at the scale of the National Soil Inventory to 
explore these effeets. 

The fact that the losses appear to be happening across both 
countries irrespective of land use suggests a link to climate change. 
Over the survey period, the mean temperature across England and 
Wales increased by about 0.5 °C and there were also changes in 
rainfall distribution 18 . Climate change will affect soil carbón turnover 
through various processes. Increases in temperature will tend to 
increase rates of organic matter decomposition by soil microbes, 
although the magnitude of this effect and differences between soils 
are uncertain 9 . The effeets of temperature will interact in a compli- 
cated way with changes in soil moisture brought about by changing 
rainfall and evapo-transpiration patterns, and changes in atmos- 
pheric C0 2 and nitrogen deposition. In freely drained soils, warmer 
drier conditions may retard decomposition of organic matter if lack 
of moisture limits soil microbes. But in wet anoxic soils, increased 


Organic carbón content (g kg -1 ) 



Figure 3 Rates of change in soil organic carbón content, grouped by 
original carbón contents and indicated land uses, a, Arable/rotational 
grass; b, permanent grass; c, non-agricultural; d, all. The ranges of carbón 
content (g kg 1 ) are, from left to right in each panel, 0-20,20-30,30-50,50- 


100, 100-200, 200-300 and >300. Bar heights indícate mean valúes; error 
bars indícate 95% confidence intervals. The derivation of equation (1) is 
described in the text. 
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Table 1 | Estimated rates of change for all original sites 


Range of 
original C org 
(g kg -1 ) 

No. sites 
originally 
sampled 

Mean 

original C org 
(g kg -1 ) 

Total a 


Total C org 0-15 cr 

ti depth* 

Rate of change 
in C org f 
(g kg 1 yr -1 ) 

Total change 
in 0-15 cm:j: 
(Tgyr -1 ) 

(km 2 ) 

(%) 

(Tg) 

(%) 

0-20 

1,061 

13.9 

26,525 

18.7 

66.4 

7.7 

0.34 

1.68 

20-30 

1,158 

24.5 

28,950 

20.5 

111.9 

12.9 

0.13 

0.66 

30-50 

1,607 

38.5 

40,175 

28.4 

214.8 

24.9 

-0.10 

-0.65 

50-100 

1,140 

66.8 

28,500 

20.1 

220.5 

25.5 

-0.68 

-2.11 

100-200 

313 

137.6 

7,825 

5.5 

92.5 

10.7 

-2.18 

-1.32 

200-300 

95 

244.5 

2,375 

1.7 

36.5 

4.2 

-4.00 

-0.59 

>300 

288 

439.7 

7,200 

5.1 

121.7 

14.1 

-7.37 

-2.10 

All 

5,662 

66.1 

141,550 

100 

864.3 

100 

-0.64 

-4.44 


* Original C org x bulk density (kgdm -3 ) x depth (dm) x area (km 2 ) x 10 —4 , where bulk density = 1.3 — (0.275 ln(C org /10)) (ref. 25). 
i'For sites not resampled, rate of change was calculated from original C org using equation (1). 

Rate of change x bulk density x area. 


oxygenation of the soil as evaporation increases the depth to the 
water table will greatly reinforce the increase in decomposition with 
increased temperature. Such effects may in part explain the increase 
in the relative rate of loss with soil carbón content because more- 
organic soils tend to be wetter. Also consistent with this, more- 
organic soils necessarily contain a greater proportion of slowly 
decaying organic matter, and the rate of turnover of this material 
appears to be more sensitive to temperature changes than that of 
more-labile organic matter 8 . However the magnitudes of these effects 
are unknown. 

On the basis of atmospheric observations, net carbón absorption 
by terrestrial systems in the Northern Hemisphere has increased in 
recent decades 14 . The magnitude of the carbón sink in different 
regions and the contributions of different processes are highly 
uncertain, but the main sinks are thought to be changes in land 
use and increased forest growth with increased atmospheric C0 2 and 
nitrogen deposition. Our findings show that losses of soil carbón in 
the UK, and by inference in other températe regions, are likely to have 
been offsetting absorption by terrestrial sinks, greatly adding to the 
uncertainty of future trends. 

METHODS 

Original sampling. The 5-km sampling grid was offset 1 km north and 1 km east 
of the origin of the National Grid to avoid sampling at the edges of published 
maps. Urban areas and water bodies were avoided, but otherwise all soils were 
sampled. At each site, 25 soil cores were taken at 4-m intervals in a 20 m X 20 m 
square centred on the grid intersection 19 . The cores were taken to a máximum 
depth of 15 cm from the surface using a 2.5-cm diameter auger, with the soil 
surface taken as the zero of measurement, and excluding litter layers and living 
vegetation. The cores were bulked in the field, giving a total of approximately 
1 kg of moist material for each site. The samples were air-dried at temperatures 
not exceeding 30 °C, and each then divided into three equal portions. One 
portion was stored in a polythene bag inside a cardboard box at ambient 
temperature and humidity without further treatment. The other portions were 
crushed to <2 mm, and subsamples analysed for organic carbón as below. 
Resampling. The minimum proportions of the original sites to be resampled so 
as to detect changes in C org with 95% confidence were calculated using the 
relation « 2/^1 — 1/[1 + ni(d/z a s) 2 ], where n x and n 2 are the sizesof the original 
and resampled populations, a is the probability that the change in mean C org is 
greater than some small valué d,z a is the probability of the standardized normal 
distribution of C org being less than ce, and s is the standard deviation of the 
original population. The accuracy of the C org measurement in the laboratory (see 
below) was ± 1 g kg - and henee the minimum valué of d is 2 g kg - *. This valué 
was used for the arable/rotational grass and permanent grass sites. For the non- 
agricultural sites, the standard deviation of C org was greater, particularly for 
more-organic soils. Accordingly, for these sites sampling was designed with 
d= 2 gkg -1 for soils with C org < 150 gkg -1 , d= 5 gkg -1 for C org — 150— 
300 gkg -1 and d= 10gkg -1 for C org > 300gkg -1 . Sites were selected at 
random from the original sites within these groupings in proportion to the 
original regional sampling density, and sampled exactly as originally. 

To test the accuracy with which the sites could be relocated, six surveyors were 
instructed to revisit 10 sites with widely differing characteristics foliowing the 
originally recorded site descriptions, their positions recorded with a global 
positioning system, and the distances between them subsequently measured. 


This showed the accuracy of relocation was better than 20 m on enclosed land 
and better than 50 m on open land, and valúes of C org at 0,20 and 50 m from the 
target measured using the original methods were not significantly different from 
each other (Supplementary Fig. 1). Hand-written field records from the original 
sampling were examined to see if outliers in the data (based on the log 
normalized changes in C org ) had any common features. No artefaets were 
found that would have led us to believe the outliers were not true representations 
of the changes. 

Organic carbón analysis. Exactly the same methods were used for the two 
samplings. Soils with C org < approximately 150gkg -1 were analysed by a 
modified Walkley-Black method 20 . The small subsamples this uses are unrepre- 
sentative of highly organic soils, so soils with C org > approximately 150 gkg -1 
were analysed by loss on ignition (LOI) 21 , converted to C org by C org = 0.5 X LOI. 
To check that the change of method introduced no artefaets, we applied both to 
95 soils with C org = 20-200 gkg -1 and obtained good agreement and no 
systematic deviation (Supplementary Fig. 2). To check for differences in 
analytical precisión between the samplings, we reanalysed stored samples from 
10 % of the original sites and obtained good agreement with the original valúes 
across the full range of C org (Supplementary Fig. 3). 

Derivation of equation (1). Equation (1) is a regression, but was not fitted by 
ordinary least squares because the original systematic sampling preeludes 
the assumption that the residuals are independent random variables 22 . 
Equation (1) was fitted as a mixed model with a spatially-dependent random 
effect and a white noise term such that, for any pair of locations x, and x-, the 
expected squared difference of the residuals, 77 (x¿) and 77 (x,), is 
E[{t 7 (x¿) — 77 (Xj)} 2 ] = 2{co + Ci(l — exp[—|x, — x ; |/a])}, where c 0 and C\ are 
the variances of the white noise and spatially dependent components, respect- 
ively, and a is a distance parameter. The model was fitted by residual máximum 
likelihood 23 . 
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